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ABSTRACT 

This analysis involves an in-depth look into the OBC attitude 
determination algorithm. A review of TRW error analysis and 
necessary ground simulations to understand the onboard attitude 
determination process are performed. In addition, a plan is 
generated for the in-flight calibration and validation of OBC 
computed attitudes . Pre— mission expected accuracies are summarized 
and sensitivity of onboard algorithms to sensor anomalies and filter 
tuning parameters are addressed. 

1 . 0 INTRODUCTION 

The Gamma Ray Observatory (GRO) (see Reference 10) is a three axis 
stabilized spacecraft scheduled to be launched into a 350-450 Km 
orbit 1990 by the Space Transportation System (STS) . The GRO 
science instruments study gamma ray sources between 0.1 to 30000 
mega-electron-volts (MeV) before they are absorbed by the Earth's 
atmosphere. The spacecraft is designed to stay inertially pointed, 
using reaction wheel control, for two weeks at a time before 
maneuvering to the next gamma ray target. 

GRO has an onboard attitude determination accuracy requirement of 
86.4 arcseconds per axis (3 sigma) during the normal science 
observation mode. This accuracy is accomplished by the use of two 
Fixed Head Star Trackers (FHSTs) and an Inertial Reference Unit 
(IRU) . Both of these attitude sensors have been used on the Solar 
Maximum Mission (SMM) , LANDSAT 4, and LANDSAT 5 spacecraft. As a 
backup, the Fine Sun Sensor (FSS) can take the place of a FHST with 
the resultant attitude accuracy of 167.5 arcseconds/axis (3 sigma). 
In both cases, the attitude is propagated using the IRU data and 
updated after a FHST or FSS measurement by using an extended Kalman 
Filter. 


2.0 GRO ONBOARD ATTITUDE ESTIMATION (Reference 1) 
Time Propagation 


In GRO, the attitude computations are contained in two modules: 
kinematic integration module and attitude estimation module. The 
kinematic integration routine uses the previous cycle OBC quaternion 
and the current gyro output to update the OBC quaternion. The 
kinematic equation for updating a quaternion is (Reference 2): 


^n+l) = <cos(wT/2) I + l/w sin(wT/2) 
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where w = ( + w~ ) 


1/2 


^ y 2 
T = time interval 
I = 4x4 identity matrix. 


Since the gyro output consists of three angles 0 X , 0y, 0 Z the 
following substitutions can be made: 


®x = w x T ' e y = w y T ' e z = w z T 

Equation 2.1 then becomes: 


q(t n+ i) = j cos (0/2 ) I + l/0sin(0/2) 
where 0 = ( © x + ©y + ©z ) 
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Every 32.768 seconds, the attitude estimation routine (ATTEST) 
generates roll, pitch, and yaw errors. These errors are fed into 
the kinematic integration routine in place of the normal gyro data 
that is used between 32.768 second updates. 


The attitude estimation routine (ATTEST) itself consist of an 
extended Kalman filter (KF) . Reference 1 contains an outline of 
ATTEST. The KF is implemented in two steps. First, the propagation 
of the internal statistics based on the Dynamics Model and second, 
updating the state vector based on the Observation Model , the 
measurements, and the internal statistics. ATTEST alternates 
between the two sensors (FHST/FHST or FHST/FSS) being used for 
attitude estimation every 32.768 seconds. 


Dynamics Model 

The gyro rate measurement is assumed to have the following form: 
© = w + b Q + b - n v 

fc = n u 

where, 0 - gyro rate measurement 
w - true spacecraft rate 
b Q - gyro bias error 
b° - gyro random walk error 

n - float torque noise (gaussian white noise) 

n^- float torque derivative noise (gaussian white noise) 

The gyro drift error, e, is defined as follows: 

e = w - © 

It then becomes the following equation: 
e = -b Q - b + n v 
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The gyro bias is assumed to be known and can be taken out of the 

above equation. Therefore: 

e = -b + Bv 

= n u 

The attitude error, is computed as follows: 
if + w x ]|J = e 

However since w is neglible, the dynamic model is reduced to the 
following form: 

$ = -b + IV 

fe = 

If these two equations are put into a linear state space formulation 
equations (2.3) and (2.4) are derived: 

X(t) = F X(t) + W(t) (2.3) 

°3x3 “ X 3x3 $ — v ( 3x1) 

+ (2.4) 

_°3x3 °3x3 _ u (3x1) _ 

where, lj[ - attitude error 

b - gyro random walk error 

n v - float torque noise (Gaussian) 

n u - float torque derivative noise (Gaussian) 

The state equation is discretized to the following form: 

X(t k ) - 0k X^k-l) + H(tk> 

where p k = e AT k and T k = t k - t k-1 . 

The two characteristics of W(t) are the mean: 

E [ W(t) ] = 0 

And the covariance: 

m ®3x3 n 

E [ W(t) W T (t') ] = JX ^ /(t-t') (2.5) 

_°3x3 ISu n U _ 

where "T" denotes the transpose. Note that the off diagonal 
elements in (2.5) are zero since it is assumed that there is no 
correlation between n u and n y . 

The Spectral Density Matrix is defined as follows: 

Q(t) = E [ W(t) W T (t) ] (2.6) 


$ 

X(t) = . 

b 
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(2.7) 


Thus the covariance is given as 

Q (t) <f(t T t') 

The discrete Dynamics Noise Covariance matrix, Q k , is obtained using 
the state transition matrix, p k , and the Spectral Density matrix, 
Q(t), in the following manner: 

Q k =^J9(t k ,t') Q(t') 0 T (t k ,t') dt' (2.8) 

fc k-l 

Once Q k is computed, it is used to propagate the state covariance 
matrix as follows: 

p k <“> = ?k p k-l< + > ?k + Q k (2,9) 

where P v (-) is the Propagated Covariance Matrix at time k and, 
p k-l is the u Pd ated Covariance Matrix at time k-1. 

Observation Model 


In the GRO Flight Software, the FHST measurements are used to 
create an observed star unit vector, OS, in the Sensor Coordinate 
frame. The identified star position in the star catalog is used to 
create an expected or computed unit star vector, CS, in the Sensor 
Coordinate Frame. We then define 


Z k (i) = OS k (i) - CS k (i) for i=x and y 
where Z k is the measurement residual. 

From this definition of Z k , H k is shown to be 


( 2 . 10 ) 


H k = 


(X X S k ) 
_( X x S k )' 


J lx3 
J lx3 _ 


( 2 . 11 ) 


where S v is the observed star vector in the spacecraft body frame, 
X is the X axis of the FHST in the spacecraft body frame, 
and Y is the Y axis of the FHST in the spacecraft body frame. 


In the observation model 

Zk = H k X + V k (2 ‘ 12 

where Z k is the observation defined in (2.9), and V k is the sensor 
noise (Gaussian) . 

The sensor noise characteristics are the following: 
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E[ v k ] = o 

R k = EC v k v^] 



(2.13) 

(2.14) 


It ^ i ® fur ^her assumed that the initial state vector X_ is Gaussian 
and Xq , W ,and V k are independant of each other. Since all are 
assumed zero mean and Gaussian, this is equivalent to assuming they 
are uncorrelated with each other. 9 y 


FSS Model (Reference 10) 

As FHST ' FSS model uses an observed Sun position, OS 

and a computed Sun position, CS, to compute measurement residual!, 
7, . as fnl 1 nwa • ' 


Z^fi) - OS k (i) - cs^fi) for i=x and y. 


The measurement equation is the same used for the FHST (2.12). For 
the FSS, the is shown to be: 


H k = 


( X M p X S k ) T o lx3 
_( Y M P x Sk ) T 0 lx3. 


(2.15) 


where S k is the computed sun vector. X MP and Y Mp are as follows 
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where X F is the FSS X-coordinate axis in the spacecraft frame, 

Yp is the FSS Y-coordinate axis in the spacecraft frame, 

Z F is the FSS Z-coordinate axis in the spacecraft frame" 

XP and YP are the FSS expected measurements. 


The FSS noise characteristics are the same as those for the FHST. 
Update Algorithms 


T ^fw S ^ te . Vector is u P dated b Y processing the following equation 

oh«^- lnpUt % Pk( " ) (2>9) ' H X < 2 * u and 2 ' 15 )' R k (2.14), and the 
observation vector Z k (2.12): K 


&<♦: k LW’ + 

X k <*> - **<-> . K>J 2j * H k 2k <-) ) 


-1 


(2.16) 

(2.17) 

(2.18) 


where P k (+) is the updated Covariance Matrix. 
K k is the Kalman Gain Matrix. 

X k (+) is the updated State Vector. 


i 
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The GRO Flight Software employees a scaler implementation method 
which requires the sequence of equations (2, 16-2.18) to be executed 
twice. In the first pass the following substitutions are made for 
the FHSTs (similarly for the FSS) : 


H k = H k,l = t ( X x s k ) T o lx3 ] 

R k = R k’l = R 11 


The resulting Kalman Gain Matrix K k , = K k is used to 
covariance matrix (2.17) where P k , = Pv and the update 
£ (_) = o. The equations are listed below: 


p k'i( + ) = i 1 - K,i H k,i3 p k(-> 

x k'l (+) = K k,l z k, 1 


p k(-> H k,l 


<k,l 


update the 
(2.18) where 


] 


where Z k ^ i s the X component of Z k . 

In the second pass, there are the following substitutions: 


Hir = H 

R £ = R 


K v = K 


k, 2 
k, 2 


= [ ( X X S k ) 


= R 
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?k(") = 


k, 1 


(+) 


x k<-) = ^k;i( + ) 


°1X3 ^ 


where X k ^ ( + ) is the state vector update from the first pass. 

The final Kalman Gain Matrix K k = K k 2 is used to update the 
covariance matrix and the state vector. The equations are listed 
below: 


K k 2 = p k("> H k / t H k p 
Pkt + ) = J- I - K k H k ] P k ?-) 


P k (-) 


X k (+) = X k (-) +”K k r Z k ; 2 - H k 
where Z k 2 is the Y component of Z k . 


H k + R k^ 
X k (-) ] 


3 . 0 ADEAS 


The Attitude Determination Error Analysis System (ADEAS) was the 
attitude tool used in this analysis. It models state estimation 
using either a batch filter or a Kalman filter. The ADEAS Kalman 
filter is the same two pass filter implementation as described for 
GRO in Section 2.0. Thus, ADEAS provides a convient method for GRO 
onboard attitude error analysis. 


ADEAS models an attitude system by use of consider and solve for 
parameters. The solve-for parameters are those the onboard filter 
uses in its state vector. For GRO, the solve-for state components 
are the three attitude errors and the three gyro drift errors. The 
consider parameters are those that the onboard filter does not take 
in account such as misalignments and scale factor errors. 
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4.0 ONBOARD ATTITUDE ACCURACY 
4.1 Attitude Accuracy Requirement 

Using the above attitude estimation algorithm, the overall attitude 
?on™i n ?Re?SrenS 1 ^“ ent "* ^ Stated capabiliti es are as 


Attitude determination accuracy 
using two FHSTs 

Attitude determination accuracy 
using 1 FHST and 1 FSS 


Requirement 

(arc-sec) 

(3-sigma) 

86.4 

167.5 


Capability 

(arc-sec) 

(3-sigma) 

71.2 

143.8 


The FHST/FHST algorithm errors were the largest single contributer 
® r £° r budget at 53.3 arcseconds/axis. The FSS alignment 

FHST/FSS prror- a hnH^ ndS { aX1S ^ WaS the lar 9 est contributer to the 

attitude deteL^ d ?® fc ' ACC ? rdlng to TRW error budgets, the absolute 
attitude determination requirement is met for both cases with a 
reserve. 

ADEAS simulations were conducted to independantly verify that the 
t^TPw d atblt V de determination requirement could be met. According 
to T RW analysis, the update filter required 3 hours to converge 

3 l : Th f slmulatlon length consisted of the three hour 

convergence time plus one orbit of data. The primary error 
parameters used in these simulations are listed below: 


Input parameters 

Dynamic Noise (Reference 3) 


4 . 2459E-2 arcseconds/sec 1 / 2 per axis 
Random walk 4.4413E-5 arcseconds/sec 3 / 2 per axis 
Misalignments 

FHST # 1 32 arcseconds/axis (Table 1) 

FHST #2 32 arcseconds/axis (Table 1 ) 

FSS 97.5 arcseconds/axis (Table 2) 

Gyro 56 arcseconds/axis (Reference 4) 

Measurement Noise 

FHST #1 32.3 arcseconds (Note 1) 

FHST #2 32.3 arcseconds (Note l) 

FSS 82.1 arcseconds (Note 2) 

Gy^g Sc ale Factor Error (3000 PPM) (Reference 5) 
l n ^^ a l Attitude Error (1800 arcseconds/axis) 

Jmtial g yro Drift Error (0.5 arcseconds/second/axis) 
Systematic Calibration Errors (FHST #1 and FHST #2) 

H and V measurements - 7 arcseconds (Reference 6) 
a ttitude Stabilization Errors 0.096 degrees (Reference 3) 
Simulation Length (273.58 minutes) 

For each simulation, the algorithm uncertainties and the jitter due 
to reaction wheel disturbances were RSS'd with the resultant 
simulation attitude errors in order to compare the simulation 
results to the error budget. 
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4 . 2 FHST/FHST CASE 


This case is the nominal configuration for attitude estimation 
onboard and the most accurate. After the three hour convergence 
period, the maximum filter roll, pitch, and yaw attitude errors were 
64 11 64.74, and 60.80 arcseconds respectively. The gyro drift 

errors for the roll, pitch, and yaw axes were 3.335x10 , 

3 356xl0~^ 3.899x10”^ arcseconds/second respectively. The onboard 

attitude accuracies are well within the 86.4 arcsecond/axis 
requirement. From Figure 4.1, steady state convergence occurs 
approximately 6200 seconds into the simulation. 


4 . 3 FHST/FSS Case 

This case is used only if one FHST fails, and the resultant accuracy 
degrades considerably. After the three hour convergence Period, the 
maximum roll, pitch, and yaw attitude errors were 126.85, 120.13, 
and 80.12 arcseconds respectively. The roll, pitch, an^ yaw gyro 
drift errors were 4.046xl0 -3 , 3.379x10 and 6.376x10 
arcseconds/second respectively. The attitude accuracies are well 
within the 167.5 arcsecond/axis requirement. From Figure 4.2, the 
steady state convergence occurs at approximately 8800 seconds. 

4.4 1 FHST With Two Guide Stars Case 

A simulation of 56000 seconds was made using FHST #2 and two guide 
stars within one degree of the FHST boresight. The stars were 
measured alternately every 32 seconds. The attitude estimation 
errors were smaller for this case than for the FHST/FSS case. From 
Figure 4.3, the convergence time for this simulation was 
approximately 48000 seconds which is about eight times that of the 
FHST/FHST case and six times that of the FHST/FSS case. The longer 
convergence time is understandable from observability reasons alone. 
The maximum roll, pitch, and yaw attitude estimation errors over the 
last 8000 seconds of the simulation were 100.34, 100.03, and 64.38 
arcseconds respectively. The roll, pitch and yaw gyro drift error 
were 5.378X10 -5 , 5.375X10 -3 , and 3.382x10 J arcseconds/second 
respectively . 

5.0 SENSITIVITY ANALYSIS 

For the sensitivity analysis, the consider parameters used for the 
ADEAS simulations were increased by a factor of two. The resul ing 
ADEAS attitude errors are RSS'd with the algorithm implementation 
error and jitter due to the reaction wheels. The case designations 
were as follows: 

1 From reference 7, the FHST calibration error is 30 arcseconds 

(3-sigma) , and the noise equivalent angle (NEA) is 24 arcseconds 
(3-sigma). The NEA is reduced to 12 arcseconds by data averaging 
onboard. The resultant measurement error is the RSS of 30 and 
12 arcseconds . 

2 From reference 8, the calibration error is 79.2 arcseconds 

(3-sigma), and the noise equivalent angle is 21.6 arcseconds. 

The resultant measurement error is the RSS of 79.2 and 21.6 
arcseconds . 
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Case 

1A 

2A 

3A 

4A 

5A 

6A 

7A 

8A 


Case 

IB 

2B 

3B 

4B 

5B 

6B 

7B 

8B 

9B 

10B 


FHST / FHST Conf igurat ion 
Error Description 
Basel ine 

2 x Gyro white noise about each gyro axis 
2x Gyro random walk noise about each gyro axis 
2x Gyro scale factor error about each gyro axis 
2x Gyro misalignment about each gyro axis 
2x FHST #1 misalignment about each FHST #1 axis 
2x FHST #2 misalignment about each FHST #2 axis 
2x FHST noise 

FHST/FSS Conf igurat ion 
Error Description 
Baseline 

2x Gyro white noise about each gyro axis 

2x Gyro random walk noise about each gyro axis 

2x Gyro scale factor error about each gyro axis 

2x Gyro misalignment about each gyro axis 

2x FHST #1 misalignment about each FHST #1 axis 

2x FSS misalignment about each FSS axis 

2x FHST noise 

2x FSS noise 

2x FHST and FSS noise 


Case 

1C 

2C 

3C 

4C 

5C 

6C 

7C 


FHST with Two Guide Stars Configuration 
Error Description 
Baseline 

2x Gyro white noise about each gyro axis 
2x Gyro random walk noise about each gyro axis 
2x Gyro scale factor error about each gyro axis 
2x Gyro misalignment about each gyro axis 
2x FHST #2 misalignment about each FHST axis 
2x FHST noise 


.tude 

error results for the 

FHST/FHST case are as 

Attitude 

Determination 

Errors 

(arcseconds) 

Case 

Axis 

Roll 

Pitch 

Yaw 

1A 


64.11 

64.74 

60.80 

2A 

X 

64.24 

64.75 

60.80 

2A 

y 

64.23 

64.86 

60.80 

2A 

z 

64.11 

64.74 

60.92 

3A 

X 

65.39 

64.80 

61.22 

3A 

y 

64.18 

66.04 

61.22 

3A 

z 

64.11 

64.74 

61.74 

4A 

X 

64.11 

64.74 

60.80 

4A 

y 

64.11 

64.74 

60.80 

4A 

z 

64.11 

64.74 

60.80 

5A 

X 

64.11 

64.74 

60.80 

5A 

y 

64.11 

64.74 

60.80 

5A 

z 

64.11 

64.74 

60.80 

6A 

X 

74.80 

76.02 

61.22 

6A 

y 

64.00 

64.74 

74.35 

6A 

z 

64.12 

64.74 

60.80 
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7A 

X 

74.77 

76.05 

61.22 

7A 

y 

64.12 

64 .,73 

75.42 

7A 

z 

64.12 

64.74 

60.81 

8A 


65.39 

66.10 

60.47 

gyro drift error 

results 

for the FHST/ FHST case are as follows 

Gyro 

Drift Estimation 

Errors ( 10~ 3 

arcseconds/second) 

Case 

AXIS 

Roll 

Pitch 

Yaw 

1A 


3.335 

3.356 

3.899 

2A 

X 

3.414 

3.352 

3.899 

2A 

y 

3.332 

3.432 

3.899 

2A 

z 

3.335 

3.356 

3.982 

3A 

X 

6.228 

3.404 

4.414 

3A 

y 

3.385 

6.260 

4.414 

3A 

z 

3.335 

3.356 

6.185 

4A 

X 

3.335 

3.356 

3.899 

4A 

y 

3.335 

3.356 

3.899 

4A 

z 

3.335 

3.356 

3.899 

5A 

X 

3.335 

3.356 

3.899 

5A 

y 

3.335 

3.356 

3.899 

5A 

z 

3.335 

3.356 

3.899 

6A 

X 

3.356 

3.374 

4.414 

6A 

y 

3.165 

3.201 

6.250 

6A 

z 

3.335 

3.356 

3.899 

7A 

X 

3.356 

3.374 

4.414 

7A 

y 

3.356 

3.376 

5.774 

7A 

z 

3.335 

3.356 

3.899 

8A 


3.805 

3.859 

3.420 


increased gyro white noise and random walk noise about an axis 
primarily affects that axis as expected. There is some correlation 
between the X and Y axes but not enough to be significant. Of the 
two errors, the random walk component proves to affect the attitude 
errors the most. As expected, the random walk errors contribute the 
most to the gyro drift estimation errors (see equations 2. 4-2. 9). 
Between updates, this higher gyro drift error would degrade the 
attitude solution since the gyro data would compensated with an 
incorrect gyro drift estimate. The gyro scale factor errors and 
misalignments have no significant affect on the attitude solution 
since the spacecraft is inertially pointed and has no significant 
angular rates. As expected, the FHST misalignments have the largest 
affect on attitude accuracy. For both FHST #1 and #2, the 
misalignment of the X and Y tracker axes result in attitude 
estimation errors of over 70 arcseconds. Both FHSTs had a maximum 
attitude estimation error of 76 arcseconds when their Y-axis was 
misaligned. The increased FHST noise has only a small effect on the 
attitude error since there are sufficient measurements to reduce the 
scope of the error, and the system has good observability. 

The attitude error results for the FHST/FSS cases are as follows: 

Attitude Determination Errors (arcseconds) 

Case Axis Roll Pitch Yaw 


331 



IB 


127.21 

2B 

X 

127.24 

2B 

y 

127.24 

2B 

z 

127.21 

3B 

X 

128.04 

3B 

y 

127.43 

3B 

z 

127.21 

4B 

X 

127.21 

4B 

y 

127.21 

4B 

z 

127.21 

5B 

X 

127.21 

5B 

y 

127.21 

5B 

z 

127.21 

6B 

X 

149.19 

6B 

y 

127.21 

6B 

z 

127.21 

7B 

X 

208.19 

7B 

y 

127.24 

7B 

z 

127.21 

8B 


128.35 

9B 


131.64 

10B 


132.30 


121.30 

74.65 

.121.30 

74.65 

121.40 

74.65 

121.30 

74.69 

120.13 

80.12 

122.02 

74.65 

121.30 

74.91 

121.30 

74.65 

121.30 

74.65 

121.30 

74.65 

121.30 

74.65 

121.30 

74.65 

121.30 

74.65 

121.30 

74.66 

121.30 

91.91 

121.30 

74.65 

209.49 

74.76 

121.36 

74.85 

121.33 

74.65 

120.32 

94.66 

126.30 

74.93 

126.36 

74.97 


The gyro drift error results for the FHST/FSS case are as follows: 


Gyro 

Drift 

Estimation Errors 

( 10 3 

arcseconds/secc 

Case 

Axis 

Roll 

Pitch 

Yaw 

IB 


4.093 

4.187 

4.169 

2B 

X 

4.147 

4.183 

4.169 

2B 

y 

4.108 

4.262 

4.169 

2B 

z 

4.093 

4.187 

4.248 

3B 

X 

7.139 

3.399 

6.376 

3B 

y 

4.277 

6.653 

4.169 

3B 

z 

4.093 

4.187 

5.962 

4B 

X 

4.093 

4.187 

4.169 

4B 

y 

4.093 

4.187 

4.169 

4B 

z 

4.093 

4.187 

4.169 

5B 

X 

4.093 

4.187 

4.169 

5B 

y 

4.093 

4.187 

4.169 

5B 

z 

4.093 

4.187 

4.169 

6B 

X 

4.093 

4.187 

4.169 

6B 

y 

4.093 

4.187 

4.486 

6B 

z 

4.093 

4.187 

4.169 

7B 

X 

4.093 

4.187 

4.172 

7B 

y 

4.093 

4.187 

6.534 

7B 

z 

4.093 

4.187 

4.172 

8B 


5.195 

3.838 

5.767 

9B 


6.073 

6.268 

2.968 

10B 


6.343 

6.329 

3.533 


or random walk about an axis primarily affects that axis. Not only 
does a correlation of the X and Y axes exist as in the FHST/FHST 
case, but a X and Z correlation exits. Of the two gyro noises, the 
random walk error has the largest affect on the attitude errors due 
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to attitude computations between measurement updates as mentioned in 
the previous case. The gyro scale factors and misalignments have no 
significant affect on the attitude or gyro bias errors due to the 
spacecraft being inertially pointed. The FHST #2 and FSS 
misalignments are the greatest contributors to attitude errors as 
expected since they define the attitude. The FSS affects the 
attitude errors more due to the greater alignment errors as compared 
to FHST # 2. Increasing the FHST noise results in the roll error 
increasing by 1 arcsecond, the pitch error decreasing by 1 
arcsecond, and the yaw estimation error increasing by 20 arcseconds. 
The increased FHST noise to 64 arcseconds is much closer to the FSS 
measurement noise of 82 arcseconds. Thus, the FSS measurements 
would be weighted almost as much as the FHST measurements. The 
large FSS alignment uncertainties are then fed into the yaw 
estimation errors. The above maximum error occurs when the FHST is 
occulted which further degrades the yaw solution. Increasing the 
FSS noise increases the roll, pitch, and yaw estimation errors due 
to the limited memory of the filter to measurements. Thus, the 
measurement noise cannot be averaged out totally. The yaw error is 
only slightly higher than the baseline since any FSS is weighted 
much less than in the baseline while the FHST measurements are 
weighted the same. When the noise on the FSS and FHST are 
increased, the total error is due totally to the increased noise and 
the limited memory of the filter mentioned above. 

The attitude error results for the 1 FHST with 2 Guide Stars are: 


Case 

Attitude Determination Errors 
Axis Roll Pitch 

(arcseconds) 

Yaw 

1C 

2C 

X 

100.34 

100.43 

100.03 

100.07 

64.38 

64.38 

2C 

y 

100.40 

100.13 

64.38 

2C 

Z 

100.34 

100.03 

64.56 

3C 

X 

109.61 

107.65 

64.38 

3C 

y 

108.02 

109.23 

64.38 

3C 

z 

100.34 

100.03 

65.98 

4C 

X 

100.34 

100.03 

64.38 

4C 

y 

100.34 

100.03 

64.38 

4C 

z 

100.34 

100.03 

64.38 

5C 

X 

100.34 

100.03 

64.38 

5C 

y 

100.34 

100.03 

64.38 

5C 

z 

100.34 

100.03 

64.38 

6C 

X 

107.77 

107.46 

64.38 

6C 

y 

100.34 

100.03 

84.97 

6C 

z 

107.77 

107.43 

64.38 

7C 

drift error 

155.73 

results 

154.99 

for the 1 FHST 

65.19 
with 2 Guide 
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Case 

Axis 

Roll 

Pitch 

Yaw 

1C 


5.378 

5.375 

3.382 

2C 

X 

5.407 

5.393 

3.382 
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2C 

2C 

3C 

3C 

3C 

4C 

4C 

4C 

5C 

5C 

5C 

6C 

6C 

6C 

7C 


y 

z 

X 

y 

z 

X 

y 

z 

X 

y 

z 

X 

y 

z 


5.396 

5.378 

8.856 

6.660 

5.378 

5.378 

5.378 

5.378 

5.378 

5.378 

5.378 

5.378 

5.378 

5.378 

6.811 


5.400 

5.375 

6.653 

8.845 

5.375 

5.375 

5.375 

5.375 

5.375 

5.375 

5.375 

5.375 

5.375 

5.375 

6.804 


3.382 

3.503 

3.382 

3.382 

6.394 

3.382 

3.382 

3.382 

3.382 

3.382 

3.382 

3.382 

3.382 

3.382 

3.694 


As in the FHST/FHST case, an increases in the white and random walk 
noise about an axis primarly affects that axis. The correlation 
between the X and Y axes still exists. Of the two noises, the 
random walk has the largest affect on the attitude due to the 
increased gyro drift estimation error as mentioned above. The gyro 
scale factor error and misalignments as expected have no affect on 
the attitude estimation error. Also as expected, the added FHST 
alignment errors have a significant affect on the attitude 
estimation errors. The maximum attitude estimation error from the 
FHST misalignments is 107 arcseconds. The primary source of 
attitude estimation error is from the FHST noise. This results from 
the small separation of the guide stars in the FHST (2.8 degrees) . 
The attitude estimation errors due to the measurement noise is an 
arctangent relationship seen in the following diagram (Reference 9) : 



For this simulation, this error was approximately 1290 arcseconds. 
The RSS attitude estimation error for the noise simulation is 
approximately 229 arcseconds. The filter was able to improve the 
solution by 1061 arcseconds. The filter's memory for the 
measurements was not long enough to average out the measurement 
noise further. 

Of the three cases simulated, the FHST/FHST case is by far the most 
accurate as expected. The choice for the backup case is not as easy 
to choice. The disadvange of the FHST/FSS case is the FSS alignment 
errors. Unless the FSS alignment is updated often, alignment errors 
due to thermal effects could degrade the attitude accuracy 
significantly. From an alignment point of view, the 1 FHST case is 
the most preferable due to the FHST being mounted on a fairly stable 
platform which is beside the gyros. From a noise point of view, the 
1 FHST case is worse than the FHST/FSS case. The accuacy could be 
improved if more than two guide stars were available or the angular 
separation was larger. A possible solution to the problem is to use 
the 1 FHST case as a backup provided at least two guide stars are 
available and the FHST noise has not increased significantly from 
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launch. The alignment errors would be less than that for the 
FHST/FSS case and more of the sky would be open for viewing due to 
the limited FOV of the FSS. Otherwise/ if the FHST noise has 
increased significantly since launch or there are not enough guide 
stars, use the FHST/FSS case. 

6.0 OBC ATTITUDE ESTIMATION CALIBRATION AND VALIDATION 
Nominal Operations 

After the spacecraft has completed inorbit checkout, there are two 
phases for validating the onboard calculated attitude: post normal 

maneuver mode phase and normal pointing mode phase. Once the normal 
maneuver is performed, the onboard filter convergence needs to be 
validated. This convergence is defined as steady state operation 
and should occur within 3 hours of filter initialization (Reference 
3). Once the filter has converged, the mission is in the normal 
pointing mode phase. According to Teledyne documentation, the 
acceleration insensitive drift rate (AIDR) peak over six hours is 
0.0006 arcseconds/second (Reference 4). Assuming the update filter 
was disabled, the attitude errors due to the AIDR alone could be 
51.84 arcseconds in 24 hours. This requires the onboard attitude 
estimation to be checked a minumum of once/day. 

Before the procedure for validating the onboard attitude estimation 
process is discussed, the error comparison limits need to be 
determined. The GRO requirement for attitude determination is an 
absolute requirement. This absolute requirement references the 
attitude to the spacecraft body. An attitude sensor alignment can 
be determined relative to an optical cube on the ground. However 
once the spacecraft is placed on orbit, this alignment is unknown 
due to launch shocks. A ground system can align the attitude 
sensors relative to a reference attitude sensor. The attitude can 
be determined relative to the reference attitude sensor. The 
resulting attitude estimation error would be a function of the 
attitude sensors and attitude determination algorithm accuracies. 
This '•relative” attitude estimation is what will be checked by the 
ground system since the alignment relative to body necessary for an 
"absolute” attitude estimate is unknown on the ground as well as 
onboard . 

The portion of the TRW error budget devoted to algorithm errors is 
53.3 arcseconds for FHST/ FHST case and 66.3 arcseconds for the 
FHST/FSS case. The FHST and FSS noise allocations in the error 
budget are 10 and 24.1 arcseconds respectively. After RSSing the 
appropriate sensor noises with the associated algorithm errors for 
each case, the following comparison accuracy limits are obtained: 


GRO ONBOARD ATTITUDE ESTIMATION ERROR COMPARISON LIMITS 

FHST/ FHST 55 . 1 arcseconds/axis 

FHST/FSS 71.2 arcseconds/axis 

Figure 6.1 


The onboard attitude estimate will be compared to the attitude 
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estimate from the Code 550 Fine Attitude Determination System 
(FADS). An outline for the comparison .procedure is as follows: 

1) Select an orbit of data. 

2) Process the sensor data through a batch least squares 
filter to obtain a ground attitude estimate. 

3) Compute errors between the ground attitude estimate 
and the OBC estimate. 

4) Compare these errors to the comparison numbers in Figure 

6 . 1 . . . 

5) If the attitude errors are less than those in Figure 6.1, 
then the OBC attitude estimation function is operating 
properly, and the validation process is complete. 

6) If the attitude errors are greater than those in Figure 
6.1, follow steps 7-9 since the the OBC may not be 
functioning properly • 

7) Check to see if the attitude estimate is diverging from 
the ground solution. If not, then the filter needs 
tuning. 

8) Check to see if any of the following attitude sensors have 
failed or if the sensor data is degraded: 

a) Gyros 

b) FHSTs 

c) FSS 

9) Check to see if the FHST(s) are tracking stars for less 
than 32.768 seconds. If this is happening, then there 
could be a star match problem. This could mean a bad 
onboard covariance (filter tuning probably required), 
failure of the FHST, or a bad uplink of guide stars. The 
update filter validates an observed star only if matches 
with one uplinked guide star. A lack of a guide star 
match or a match with more than guide star causes the 
update filter to send a break track command to search for 
a new star. 

At the time these tests are being performed, all update filter data 
base parameters need to be checked for the previous 24 hour period. 
Bad data base updates could easily upset onboard attitude 
estimation. 

Contingencies 

If the attitude estimation errors equal or exceed comparison limits 
and the possible problems discussed above have been eliminated, OBC 
attitude determination calibration may be required. The Onboard 
filter can be calibrated by the following methods: 

a) tuning parameter adjustment 

b) changing update frequency of covariance matrix 

c) changing measurement frequency 

Onboard filter calibration can most easily be performed by tuning 
parameter adjustment. Short of modifying the Onboard filter, the 
tuning parameters are the following: 
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a) Initial Attitude errors along the covariance 

diagonal 

b) Initial Gyro Drift error's along the covariance 

diagonal 

c) FHST noise variance 

d) FSS noise variance 

e) Gyro white noise estimate. 

f) Gyro random walk estimate. 

The initial attitude and gyro drift errors in the state covariance 
matrix occupy the first six diagonal elements. These are data base 
constants that are used to initialize the state covariance matrix 
when the filter is initialized. The more accurate these numbers 
are, the faster the filter will converge to the correct solution. 
Since the largest value of the upper lefthand 3x3 matrix are used 
for star selection criteria, a good estimate of the initial attitude 
and gyro drift errors will decrease the possibility of a 
misidentification of a star. 

The sensor noise variances are used in determining the weight of a 
particular measurement. This can be seen in the GRO Kalman gain 
equation: 

K Jc,i = p k,i(-> H k,i / r H k,i p k,i(-> H k,i + R k,i ] (6.1) 

where R k * is the sensor noise for a particular measurement. If 
R k. i increased, then the Kalman gain K k will be decreased and 
mote emphasis will be placed on the estimate X k (-) : 

X k (+) = & k (~) + K k ( Z k - H k X k (-) ) (6.2) 

This allows the filter to place emphasis on more accurate 
measurements. In the FHST/FSS case, the FHST measurements would be 
given more confidence by the filter than the FSS measurements since 
FHSTs are more accurate. The estimated dynamic noise for the filter 
(white noise and random walk) come into play with the propagation of 
the state covariance matrix. The spectral density matrix is defined 
as follows: 


Q(t) = E [ W(t) W T (t) ] (6.3) 

where W(t) is white noise vector in the state equation (2.3). Q(t) 
is used to form the discrete dynamics noise covariance matrix as 
follows: 

tv 


Q k 


fc k-l 


t') Q(t') 0 T (t k ,t') dt' 


(6.4) 


which is used in the propagation of the covariance matrix as follows 
p k(-> = *k p k-l( + ) 0k 


+ Qk 


(6.5) 



Calibration j 

Method 1 1 

Table 1 

Parameters | < 

FDF 1 

Control? J 

Comments j 

TUNING PARAMETERS 




FHST Noise Var. 

Table 37 

YES 1 

Jseful if FHST data is 
noisy or to add weight 
to the FSS measurements 

FSS Noise Var. 

Table 42 

NO 1 

Useful if FSS data is 
noisy or to add weight 
to the FHST measurements 

Initial attitude 
errors 

Table 42 

NO 

Useful for the 
convergence of the filter 

Initial gyro 
drift errors 

Table 42 

NO 

Useful for the 
convergence of the filter 

White noise 

Table 36 

YES 

Needed if gyro white 
noise increases or if 
less emphasis on the 
dynamic model is desired 

Random walk 

Table 36 

YES 

Needed if gyro random 
walk noise increases or 
if less emphasis on the 
dynamic model is desired 

Related Parameters 




Estimated attitude 
errors 

Table 42 

NO 

Useful for small onboard 
attitude quaternion 
corrections . 

Estimated gyro 
drift errors 

Table 42 

NO 

Needed to insure the 
onboard attitude is 
propagated correctly. 

Attitude quaternion 

Table 59 

NO 

Needed if onboard 
quaternion is bad 

OTHER METHODS 




Covariance Update 
Frequency Change 

N/A 

N/A 

Requires code changes 
in ATTEST and probably 
the EXEC routine. Not 
an easy modification. 

Measurement Update 
Frequency Change 

i N/A 

N/A 

Requires changes in the 
EXEC routine. Probable 
OBC loading problems. 


Table 6.1 Update Filter Calibration Method Summary 
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®J e . F ro 2 a ? atlon of the covariance matrix is used in the computation 
°£ as 0811 be see n in equation (6.1) as well as the 

f '" ar identification process mentioned above. Since the Kalman qain 
is used to determine the updated state vector (6.2), P v (-) has 
potential to undermine the state estimation process ir poorly 
computed. 1 

OTie second possible calibration method is changing the update 
frequency of the state covariance matrix (6.5)! Currently, the 

updates P k(“) every 32.768 seconds before the star 
identification process occurs. Increasing the update frequency 
could assist in state estimation as well as star identification. 

up< ? ate is embedded in ATTEST which makes 
this possibility difficult to implement due to the software mods 
needed . 


f ast calibration method mentioned deals with measurement 
updates. Stellar updates very accurately pin down attitude errors 
and gyro drift errors. Gyro drift errors affect the attitude 
between measurement updates. During periods of no measurement 
updates, the onboard quaternion would be updated in the kinematic 
integration routine which uses slowly degrading gyro data. The gyro 
data used requires an accuate estimate of the gyro drift error to 
compensate the incremental angles. Without accurate gyro drift 
compensation, the onboard attitude quaternion would slowly diverge 
from the proper attitude. Mow if measurement updates were made more 
frequently, the attitude would be compensated before it could 
degrade significantly. This modification would be easier to perform 
than changing the propagation frequency of the state covariance 

®* nce little coding changes would be required. The resulting 
OB< ? . l° adlI ?g would need to he studied to determine whether this 
modification is viable. 


7 . 0 CONCLUSIONS 

F °£.k°£ h the FHST/FHST and FHST/FSS configurations, the GRO onboard 
attitude determination accuracies can be met with significant 
margins (assuming nominal on-orbit conditions) . For the FHST/FHST 
configuration, the roll, pitch, and yaw predicted attitude 
estimation errors are 64.11, 64.74, and 60.80 arcseconds 
respectively.. The requirement is 86 arcseconds/axis. For the 
FHST/FSS configuration, the roll, pitch, and yaw predicted attitude 
estimation errors are 127.21, 121.30, and 74.65 arcseconds 
respectively. The requirement is 167.5 arcseconds/axis. For the 1 
FHST with 2 guide star case, the attitude estimation errors were 
100.34, 100.03, and 64.38 arcseconds respectively. 

As expected for the FHST/FHST and FHST/FSS cases, the attitude 
estimation accuracies are most sensitive to FHST and FSS alignment 
errors. The gyro drift errors are most sensitive to random walk 
errors in both cases, but alignment errors provide the second 
largest component of gyro drift error. These results emphasize the 
importance of proper attitude sensor alignment determination. 

The single FHST case demonstrated that onboard attitude estimation 
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rivalling that of the FHST/FSS combination is possible. However, 
this choice is highly sensitive to measurement noise, if the 
measurement noise is nominal and two guide stars are available it 
iS recoLended to use this case over the FHST/FSS case due to high 
FSS alignment errors and the FSS FOV limitation on target attitudes. 

A orocedure was outlined for GRO OBC attitude estimation validation. 
It P was pointed out that the accuracy check is a relative check and 
not- an absolute check. The attitude error comparison for the 
FHST/FHST case is 55.1 arcseconds/axis and 71.2 arcseconds/axis for 
the FHST/FSS case. 

When and if problems with the onboard attitude estimation process 
are detected, three onboard filter calibration techniques are 
available* filter tuning, increased state covariance matrix update 
frequency, and increased ATTEST frequency. Of the three techniques, 
filter tuning parameter adjustment is the easiest. In extreme 
situations where filter turning is not sufficient, an ^creased 
measurement frequency for ATTEST would require the least software 
modifications. OBC loading would need to be studied though since 
ATTEST is the largest function in the Attitude Control and 
Determination (ACAD) portion of the OBC flight software. 
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